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ABSTRACT 

Two different simple measurements of galaxy star formation rate with differ- 
ent timescales are compared empirically on 156, 395 fiber spectra of galaxies with 
r < 17.77 mag taken from the Sloan Digital Sky Survey in the redshift range 
0.05 < z < 0.20: a ratio A/K found by fitting a linear sum of an average old 
stellar poplulation spectrum (K) and average A-star spectrum (A) to the galaxy 
spectrum, and the equivalent width (EW) of the Ha emission line. The two mea- 
sures are strongly correlated, but there is a small clearly separated population of 
outliers from the median correlation that display excess A/K relative to HaEW. 
These "K+A" (or "E+A" ) galaxies must have dramatically decreased their star- 
formation rates over the last ~ 1 Gyr. The K+A luminosity distribution is very 
similar to that of the total galaxy population. The K+A population appears to 
be bulge-dominated, but bluer and higher surface-brightness than normal bulge- 
dominated galaxies; it appears that K+A galaxies will fade with time into normal 
bulge-dominated galaxies. The inferred rate density for K+A galaxy formation 
is ~ lQ~^h^ Mpc^^ Gyr^^ at redshift z ~ 0.1. These events are taking place in 
the field; K+A galaxies don't primarily lie in the high-density environments or 
clusters typical of bulge-dominated populations. 
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Subject headings: galaxies: clusters: general — galaxies: evolution — galaxies: 
fundamental parameters — galaxies: statistics — galaxies: stellar content — 
stars: formation 

1. Introduction 

Cosmic star formation activity appears to be coming to an end. All indicators of star 
formation show that the cosmic mean comoving density of star formation has been declining 
since redshifts near unity (e.g., Hogg 2001, and references therein). The bulk of stars in the 
local Universe are found in old stellar populations (Fukugita et al. 1998; Hogg et al. 2002). 
One consequence of this global decline is that galaxies in the local Universe ought to be 
shutting down their individual star formation activities. It is not understood whether the 
changes in star formation in individual galaxies are expected to be gradual or dramatic. 

Our extremely local neighborhood in the disk of the Milky Way has been forming stars 
at a relatively constant rate for the last few billion years (Barry 1988; Pardi & Ferrini 1994; 
Prantzos & Silk 1998; Rocha-Pinto et al. 2000; Gizis et al. 2002), but it appears that it 
can't continue for much longer. When the total available gas content is compared to the star 
formation rate, it appears that the Milky Way will run out of gas in the next few Gyr (e.g., 
Mooney & Solomon 1988; Blitz 1996). This is also true for the majority of local disk galaxies 
(Kennicutt ct al. 1994). In detail these calculations do not account for gas recycling, gas 
flux from plasma resevoirs, or new gas accretion, but the trend they imply for cosmic star 
formation activity is consistent with that found in cosmological observations. 

Here we begin a comparison of different optical star-formation rate measurements in 
a large sample of galaxies with spectroscopic observations in the Sloan Digital Sky Survey 
(SDSS; e.g., York et al. 2000). One star-formation rate measurement makes use of the 
strength of the Ha emission line, which is emitted in ionization regions around O and B 
stars, with hfetimes of ~ 10^ yr. Another uses the fraction of the galaxy light emitted by 
A stars, with hfetimes of ~ 10^ yr. These star-formation indicators are sensitive to stars of 
different masses and lifetimes. 

Star-formation indicators tracking stars with different lifetimes can be compared to flag 
galaxies with strong star-formation rate derivatives. In particular, comparison of our Ha and 
A-star measurements locates galaxies that have changed their star-formation rates within 
the last ~ 10^ yr. If star formation is triggered by mergers or interactions, such galaxies 
might be found to be undergoing or recovering from an interaction, and their statistics 
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might provide a measure of the merger rate. If star-formation efficiency is a strong function 
of local environment, then galaxies with negative star-formation derivatives might be tracers 
of the environmental conditions at the transition between those that are conducive to star- 
formation and those that aren't. 

One population of galaxies with strong negative star-formation rate derivatives is already 
known: "K+A" or "E+A" galaxies (Dressier & Gunn 1983; Zabludoff et al. 1996; Goto 
et al. 2003). These are galaxies combining the stellar absorption lines of A stars with 
those of old stellar populations, but, at the same time, showing little sign of current star- 
formation. K+A galaxies arc exceedingly rare — they comprise between 10^^ and ~ 10^'^ of 
the galaxy population (depending on definition) — but they are expected to provide important 
information about the conditions under which stars form and galaxies evolve. 

For example, if it is true that elliptical or bulge-dominated galaxies are formed by 
the mergers of spiral or disk-dominated galaxies, then the central concentration of stars, the 
alpha-enhanced chemical abundances (e.g., Worthey 1998), and the lack of cold gas, dust and 
star formation (e.g., Roberts & Haynes 1994) all require large, brief, centrally concentrated 
starbursts during or immediately following disk-disk mergers. Note that the abundance 
ratios provide the evidence that the starbursts must be brief. These merger events therefore 
likely go through a K+A phase for ~ 10^ yr; in this scenario, it should be possible to use 
the abundances of K+A galaxies to constrain merger rates. Indeed, it has already been 
shown that the A stars tend to be centrally concentrated in K+A galaxies and that some 
have "kincmatically hot" stellar populations like bulge-dominated galaxies (Norton et al. 
2001). On the other hand, in some K+A galaxies, there is residual gas, implying that the 
star formation event does not always exhaust the entire gas reservoir (Chang et al. 2001). 

Although any scenario in which bulge-dominated galaxies are created by mergers virtu- 
ally requires the remnants to be K+A events, an apparent excess of A stars over O and B 
stars can be produced in heavily dust-enshrouded star-formation regions from which young 
stars escape on Gyr timescales, or in mergers in which an old galaxy captures a smaller, 
younger galaxy, shutting off star formation in the process. 

In what follows, a cosmological world model with (JImj^a) = (0.3,0.7) is adopted, 
and the Hubble constant is parameterized Hq — 100 h kms~^Mpc~^, for the purposes of 
calculating distances and volumes (e.g., Hogg 1999). 
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2. Data 

The SDSS is taking ugriz CCD imaging of lO'' dcg^ of the Northern Galactic sky, 
and, from that imaging, selecting 10^ targets for spectroscopy, most of them galaxies with 
r < 17.77 mag (e.g., Gunn et al. 1998; York et al. 2000; Stoughton et al. 2002). 

All the data processing, including astrometry (Pier et al. 2003), source identification, 
deblending and photometry (Lupton et al. 2001), calibration (Fukugita et al. 1996; Smith 
et al. 2002), spectroscopic target selection (Eisenstein et al. 2001; Strauss et al. 2002; Richards 
et al. 2002), and spectroscopic fiber placement (Blanton et al. 2003) are performed with 
automated SDSS software. 

Every spectral "plate" of fiber positions includes several faint (15.5 to 18.5 mag) F stars. 
The spectra are calibrated with these F-star spectra; ie, they are multiplied by the function 
of wavelength that makes the F-star spectra match F-star spectrophotometry. Although this 
calibration procedure produces consistent calibration at the 10-percent level (as measured 
by comparisons of old stellar populations at various redshifts, Eisenstein et al. 2003), it does 
not deal carefully with the diversity of F-star spectra, or the possibihty that some of the F 
stars, even at inferred distances > 2 kpc, may be inside some part of the Galactic extinction. 
Inasmuch as the calibration stars are beyond the Galactic extinction, however, the spectra 
are effectively extinction-corrected. 

Redshifts are measured on the reduced spectra by an automated system, which models 
each galaxy spectrum as a linear combination of stellar eigenspectra (Schlegel, in prepara- 
tion). The eigenspectra are constructed over the rest-frame wavelength range 4100 < A < 
6800 A from a high- resolution spectroscopic library (Prugniel & Soubiran 2001). The cen- 
tral velocity dispersion is determined by fitting the detailed spectral shape as a velocity- 
smoothed sum of stellar spectra (Schlegel & Finkbeiner, in preparation). 

It must be emphasized that the 3 arcsec diameter spectroscopic fibers of the SDSS 
spectrographs do not obtain all of each galaxy's light because at redshifts of 0.05 < z < 0.2 
they reprsent apertures of between 2 and 7 kpc diameter. For this reason, in what 
follows, SDSS imaging data rather than spectroscopy are used to infer the global properties 
of the galaxies. 

Galaxy luminosities and colors are computed in fixed bandpasses, using Galactic extinc- 
tion corrections (Schlegel et al. 1998) and K corrections (computed with kcorrect vl_ll; 
Blanton et al. 2003). They are K corrected not to the redshift z — observed bandpasses but 
to bluer bandpasses ^-^g, °'^r and "made" by shifting the SDSS g, r, and i bandpasses to 
shorter wavelengths by a factor of 1.1 (c.f., Blanton et al. 2003, a). This means that galaxies 
at redshift z — 0.1 (typical of the sample used here) have trivial K corrections. 
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To the azimuthally averaged radial profile of every galaxy in the observed-frame i band, 
a seeing-convolved Sersic model is fit, as described elsewhere (Blanton et al. 2003a). The 
Sersic model has surface brightness / related to angular radius r by / oc exp[— (r/ro)*^^/"-*], 
so the parameter n (Sersic index) is a measure of radial concentration (seeing-corrected). At 
n = 1 the profile is exponential, and at n = 4 the profile is de Vaucouleurs. 

To every best-fit Sersic profile, the Petrosian (1976) photometry technique is applied, 
with the same parameters as used in the SDSS survey. This supplies seeing-corrected Pet- 
rosian magnitudes and radii. A A'-corrcctcd surface-brightness fio.i^ in the "'^i band is com- 
puted by dividing half the i^-corrected Petrosian light by the area of the Petrosian half-light 
circle. 

For the purposes of computing large-scale structure statistics, we have assembled a 
complete subsample of SDSS galaxies known as the NYU LSS samplel2. This subsample is 
described elsewhere (e.g., Blanton et al. 2003a); it is selected to have a well-defined window 
function and magnitude limit. In addition, the sample of galaxies used here was selected 
to have apparent magnitude in the range 14.5 < r < 17.77 mag, redshift in the range 
0.05 < z < 0.20, and absolute magnitude Mo.ij > —24 mag. These cuts left 156,395 galaxies. 

For each galaxy, a selection volume Knax is computed, representing the total volume 
of the Universe (in Mpc^) in which the galaxy could have resided and still made it 
into the sample. The calculation of these volumes is described elsewhere (Blanton et al. 
2003b, a). For each galaxy, the quantity 1/Kiax is that galaxy's contribution to the cosmic 
number density, and the quantity L/Vmax, where L is the luminosity in some bandpass, is 
that galaxy's contribution to the luminosity density in that bandpass. 

Around each galaxy in a complete subsample, samplelO, of samplel2, there is a measure 
of the overdensity 5i inside a spherical Gaussian window of radius 1 Mpc, made by 
deprojecting the angular distribution of nearby galaxies detected in the imaging at fiuxes 
corresponding to luminosities within 1 mag of L* at the target galaxy's redshift, as described 
elsewhere (Eisenstein 2003; Hogg et al. 2003). The individual overdensity estimates are low 
in signal-to-noise, but they are unbiased in the mean when averaged over sets of galaxies. A 
galaxy in an environment with the cosmic mean density has di = 0. 

There is another measure Ss on an 8 h~~^ Mpc scale, that is a measure of the three- 
dimensional comoving number density excess around each galaxy in a sphere of 8 h^^ Mpc 
comoving radius, with no correction for peculiar velocities, as described elsewhere (Hogg 
et al. 2003; Blanton et al. 2003a). Neighbor galaxies are counted, the result is divided by 
the prediction made from the galaxy luminosity function (Blanton et al. 2003b), and unity 
is subtracted to produce the overdensity estimate. The estimates Ss are higher in signal-to- 
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noise than the estimates 5i, but they are shghtly biased because they do not account for 
pecuhar velocities. 

3. Method 

An average old stellar population spectrum (hereafter called "K" because it is domi- 
nated by K-type stars) was made by averaging the spectra of bulge-dominated galaxies with 
luminosities near L* and in typical environments, as described previously (Eisenstein et al. 
2003). An average A-star spectrum ("A") was obtained by scaling and taking the mean of 
calibrated SDSS spectra of stars with types determined by eye to be near AO. These two 
spectra were normalized by scaling them to have equal total flux within the wavelength range 
3800 < A < 5400 A. 

A linear sum of the K and A spectra was fit to every SDSS spectrum in the wavelength 
range 3800 < A < 5400 A, weighting by the inverse square uncertainty in the observed 
spectrum. A 280 kms^^ region was masked out around H/? and the [OIII] lines; this mask 
size was determined by observing that it would mask out the great majority of the emission- 
line flux in the great majority of spectra. Before fltting, the A spectrum was smoothed to 
the appropriate velocity dispersion for each spectrum. Dressier & Gunn (1983) pioneered 
this method by showing that post-starburst galaxies can be modeled well by an A dwarf 
spectrum and a KO spectrum over the range of 3400 to 5400 angstroms. They found, as we 
do, that this model simultaneously describes the broad spectral continuum and the depths 
of narrow features. 

The fitting procedure returns two numbers, the amplitudes K and A of the K and A 
spectral components in the best fit. 

The Ha line fiux is measured in a 20 A width interval centered on the line. Before the flux 
is computed, the best-flt model A-l-K spectrum is scaled to have the same flux continuum as 
the data in the vicinity of the emission line and subtracted to leave a continuum-subtracted 
line spectrum. This method fairly accurately models the Ha absorption trough in the con- 
tinuum, although in detail it leaves small negative residuals, as will be shown below. The 
flux is converted to a rest-frame equivalent width (EW) with a continuum found by taking 
the inverse-variance-weighted average of two sections of the spectrum about 150 A in size 
and on either side of the emission line. 

Some general caveats apply to our measurements. The flrst is that the K-l-A flt is not 
always good; in particular, galaxies that have undergone a very recent (< 10^ yr) starburst 
can have the blue end of their spectra dominated by O and B stars, not A stars. These 
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galaxies will be assigned high A/K, because the A template fits better than the K, even 
though the A stars are not dominating the spectra. Of course, most of these galaxies also 
show strong nebular Ha emission, and, indeed, they do contain young stars. 

We rely heavily on the flux-calibration of the spectra, because A/K is sensitive to the 
color or tilt of the spectrum. If there were large (> 20 percent) absolute calibration devi- 
ations among the galaxy spectra, we might introduce some spurious interlopers. However, 
the peak of the A/K distribution of galaxies with low H«EW at A/K = is narrow and sym- 
metrical about zero. This is a strong indicator that the spectrophotometric calibration is 
wcll-bchaved. Note that wc rely on the consistency of SDSS spectrophotometric calibration, 
but not on its absolute accuracy. The typical standard linear fitting errors in A/K are 0.03 
to 0.1, roughly consistent with the scatter around A/K = 0. 

We are looking at fiber spectra, inside the central 3 arcsec diameter of each galaxy. The 
physical size of this aperture is between 2 and 7 kpc in diameter, thereby excluding light 
from the outer parts of the galaxies. Some of our galaxies classified as K-l-A could, in prin- 
ciple, have active star formation in their outskirts. Strictly speaking, the K-l-A classification 
is a classification of the star formation indicators in the inner regions of galaxies. For this 
reason, in what follows the SDSS imaging rather than spectroscopy is used to measure global 
properties of the galaxies. 

The fitting procedure doesn't take dust exctinction into account. The A/K ratio is 
sensitive to the shape of the blue end of the spectrum, so it will be reduced in a galaxy in 
which the A stars are behind significant amounts of dust. This is also true of the HaEW 
measurements, although less so because Ha is rcdward of the K-l-A fitting region. For these 
reasons, there may be some truly post-starburst galaxies not classified as K-|-A galaxies 
because their young stars are dust-enshrouded. 

The K-l-A population defined by this project docs not contain any significant AGN 
because the selection permits so little Ha emission. There may be some truly post-starburst 
galaxies not classified as K-l-A galaxies because they contain central AGN and therefore emit 
significant Ha emission. 

4. Model spectra 

The K, A, and Ha measurements were made also on model spectra made by convolving 
PEGASE (Fioc & Rocca-Volmerange 1997) instantaneous starburst models with different 
star formation histories. 
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The original instantaneous starburst history was made with solar metallicity. Although 
nebular emission is included in the model, there was no modeling of galactic winds, infall, 
or internal extinction. For each timestep, the K and A measurements were computed with 
exactly the same least-square-fitting procedure as was performed on the data. The continuum 
around the Ha line was also computed in exactly the same way as in the data, but the fiux 
of the Ha line was taken directly from the PEGASE outputs, because PEGASE returns 
integrated flux values for the emission lines rather than flux densities. 

Three model spectra and associated evolutionary tracks were computed by convolving 
the instantaneous PEGASE burst model: A galaxy forming stars at a fixed rate for 14 Gyr; a 
galaxy forming stars at a fixed rate for 10 Gyr, shutting off, and fading quiescently for 4 Gyr; 
and a galaxy forming stars at a fixed rate for 3 Gyr, shutting off, and fading quiescently for 
11 Gyr. 

In the end, because of small differences between the old model spectra and old observed 
galaxies, the final K and A values were transformed (by a linear transformation equivalent 
to a flux-conserving shear in the K vs. A plane) so that an old model galaxy shows A/K = 0, 
just like an old real galaxy. The linear transformation affects the model outputs by less than 
15 percent. 

The measurements A/K and HaEW are plotted as a function of cosmic time for the 
three models in Figure 2. 

5. Results 

The measured A/K values for the galaxies in the sample arc plotted against the measured 
HaEW values in Figure 3. The two star-formation indicators are highly correlated, as 
expected. Overplotted on the flgure are models of galaxies with star formation histories 
consisting of constant star-formation rates followed by periods of total quiescence. 

It is worthy of note that most of the galaxies appear more quiescent (i.e., have lower 
A/K and lower HaEW) than the model (labeled "14") in which stars form at a constant rate 
over the entire lifetime of the Universe. In other words, the present-day star-formation rate 
appears much lower than the cosmic time-averaged mean. 

The models in which star formation has "shut off" drop precipitously in HaEW (which 
indicates stars with lifetimes of ~ 10^ yr) and then slowly in A/K (which indicates stars 
with lifetimes of ~ 10^ yr. The models suggest that galaxies whose star-formation activities 
have shut off will evolve into the zero-star-formation population along the zero HaEW fine. 
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Indeed, Figure 3 shows a clear "spur" of galaxies along this hne. These are the "K+A 
galaxies" . 

Figure 4 shows the ratio of HaEW to A/K for samples cut in A/K. The histograms are 
bimodal, showing a population with anomalously low HaEW. These are the K+A galaxeis. 
The figure shows that K+A galaxies form an identifiable separate population. This figure 
also suggests that the line that separates the "K+A galaxies" from the rest of the population 
has a slope of (HaEW) / (A/K) = 5 on Figure 3. This figure also shows that by defining "K+A 
galaxies" to have A/K > 0.2 we are excluding a large number of transition objects. The cut 
was put at 0.2 to avoid interlopers with true values of A/K pa 0.0 but scattered high by 
calibration or statistical errors. 

The selection of K+A galaxies for what follows is indicated by the green lines on Figure 3. 
There are 1194 K+A galaxies within these selection boundaries in sample 12. Spectra of a 
randomly selected subsample of this population is shown in Figure 5. 

The astute reader will notice that the mode HaEW is less than zero; this offset is due to 
the crudeness of the continuum model subtracted during Ho; line measurement. This offset 
does not affect any of our conclusions. 

It is interesting to compare our selection criteria with those of Zabludoff et al. (1996) and 
Goto et al. (2003). Zabludoff et al. (1996) considered galaxies in the redshift range 0.05 < z < 
0.13 with [Oil] emission EW < 2.5 A and H7, E6, and E(3 absorption EW > 5.5 A. They 
found 21 K+A galaxies out of a sample of 11113 (roughly 0.2 percent). Goto et al. (2003) also 
worked with an SDSS sample of galaxies, but in the redshift range 0.05 < z < 1.0 and with 
H(5 > 4 A. They found 3340 HDS galaxies out of a sample of 95479 (roughly 3.4 percent). 

Figure 6 shows the luminosity function of the K+A galaxies compared with the lumi- 
nosity function of the entirety of samplel2. The shapes are very similar. Interestingly, the 
shape of the K+A luminosity function appears closer to the shape of the total luminosity 
function than it does to the early- type galaxy or red galaxy luminosity functions (e.g., Blan- 
ton et al. 2003b). In detail, quantitative comparisons of the luminosity functions ought to 
involve the fading with time of the K+A galaxy light. 

Figure 7 shows the distribution of A/K for the K+A galaxies, in number-density units 
(the "ratio function"?). Overplotted are the expected distributions of A/K under the as- 
sumption of steady-state creation of K+A galaxies along evolutionary tracks like those of 
the "10" model (i.e., 10 Gyr of constant star formation followed by nothing) or the "3" model 
(i.e., 3 Gyr of constant star formation followed by nothing). Neither of these models is a 
good description, suggesting that there is probably a large diversity to the star-formation 
histories that end in a K+A phase according to the selection criteria used here. 



-10- 



Figure 8 shows the luminosity-density weighted distribution of all galaxies in the plane of 
^'^{9 ~ "f^) color and Sersic index n (seeing-corrected concentration). As discussed elsewhere 
(Blanton et al. 2003a), the general galaxy population is bimodal in this plane, showing 
a population of blue, exponential (disk-dominated) galaxies and a separate population of 
red, concentrated (bulge-dominated) galaxies. Figure 8 also shows that galaxies chosen to 
have high star-formation rates, by either A/K or HaEW, lie predominantly in the blue, 
exponential class. Galaxies chosen to have low star-formation rates lie predominantly in 
the red, concentrated class. Galaxies chosen to have high A/K but low HccEW — the K-l-A 
galaxies — have some members clearly in the blue, exponential class, and some in an outlier 
class that appears concentrated, but blue. If star formation has truly been shut off in these 
galaxies, they are expected to fade and redden and become part of the red, concentrated 
class. 

Figure 9 is similar to Figure 8 but shows the distribution of objects in the color-surface- 
brightness plane. The population of objects in this plane is not clearly bimodal, however 
the bulge-dominated galaxies form a strong peak at high surface brightnesses and red colors. 
Just like Figure 8, this shows that galaxies with high A/K or high HaEW appear to be disk- 
dominated and galaxies with low A/K or low HaEW appear to be bulge-dominated. The 
A-l-K population shows even higher surface brightnesses than the bulge-dominated galaxies, 
suggesting that as they age, the A-l-K galaxies will fade and redden into the bulge-dominated 
part of the diagram. 

Three-color images made from the SDSS imaging were inspected by eye for a randomly 
selected sample of 160 of the K-l-A galaxies in the sample. The vast majority of them 
(~ 90 percent) appear bulge dominated. The ~ 10 percent that appear to be disk dominated 
do also show significant bulges. About 30 percent either show tidal features or lie extremely 
close on the sky to nearby neighbors, indicating possible past or current interactions. 

In Table 1, the mean ovcrdensities {61) and {Sg) on the Ih^^ and 8h~^ Mpc scales are 
given for the different galaxy populations mentioned above. Error bars on these averages were 
computed by splitting the sample into three disjoint sky regions. This technique produces 
conservative error estimates because it effectively includes contributions from variations in 
calibration and large-scale structure. The numbers of galaxies used for mean ovcrdensities 
are smaller than used in the Figures because ovcrdensities have only been calculated on 
an unbiased subsample of the sample we are using. Table 1 shows that K-|-A galaxies, 
on average, live in the lower-density regions more typical of spiral galaxies, and not the 
higher- density regions typical of bulge-dominated galaxies. 
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6. Discussion 

Two star-formation indicators, HaEW and A/K, with two different timescales, ~ fO^ yr 
and ~ f 0^ yr, have been measured in a sample of f 56,395 SDSS (optically selected) galaxies in 
the redshift range 0.05 < z < 0.20. The two star-formation indicators are strongly correlated, 
as expected. Comparison with models shows that the average levels of star formation in these 
galaxies is below the average over cosmic time. This fits in with many observations of more 
distant galaxies that have implied a rapid decrease in the global comoving star formation 
rate with cosmic time to the present day (Hogg 2001, and references therein). 

A clearly separated, distinct population of "K-l-A galaxies" is found, with very weak or 
no HaEW, but significant A/K. The star-formation indicators HaEW and A/K have very 
different timescales over which, in effect, they average recent star formation rates. Galaxies 
that abruptly (ie, on timescales < 10^ yr) shut off star formation will fade rapidly in HaEW, 
but slowly in A/K. This timescale mismatch explains the presence of the distinctly visible 
K-l-A population shown in Figure 3; we expect these galaxies to indicate abruptly terminated 
star formation. 

Bulge-dominated and disk-dominated galaxies can be separated by their colors, con- 
centrations (Sersic indices), and surface brightnesses. By these measures, most of the K-|-A 
galaxies appear to be bluer, higher surface-brightness counterparts of bulge-dominated galax- 
ies. The offsets in color and surface brightness are consistent with the fading and redenning 
expected in the evolution of stellar populations. Taken on face value, this is evidence that 
K-l-A galaxies are plausibly the post-starburst progenitors of typical bulge-dominated galax- 
ies, as might be expected in a scenario in which disk-dominated galaxies, for instance, merge 
to form bulge-dominated galaxies. 

Interestingly, if that conclusion is correct, then these K-l-A galaxies are the progenitors 
of the bulge-dominated galaxies of the field, not of rich clusters, because the K-l-A galaxies 
live in the typical overdensities of spirals, not ellipticals. Past work on K-l-A galaxies has 
implied that they form a special cluster population (eg. Dressier & Gunn 1983), but this is 
presumably a consequence of the selection in those studies. 

In order to put the steady-state predictions onto Figure 7, it was necessary to multiply 
the time spent at each value of A/K, or |di/d(A/K)| by a comoving rate density to match 
the total abundance. This rate density is 1.8 or 1.0 times Mpc~^ Gyr~^ for the "3" 

or "10" model respectively. In other words, the comoving K+A galaxy creation rate density 
is on this order at redshifts z ~ 0.1. This suggests that of order 1 percent of galaxies are 
currently K-|-A, and of order 10 percent have been through a K-l-A phase since redshift unity. 

The A and K spectra used for this study will be made available with the electronically 
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published version of this article. 
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Fig. 1. — The normalized model A-star spectrum ("A") and average old stellar population 
spectrum ("K") used in the K+A fits. The vertical lines indicate the region over which 
the fit was performed. As stated in the text, the spectra have equal areas inside the fitting 
region. 
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Fig. 2. — The indicators A/K and Ha plotted against time for the three star formation 
models with constant star formation for 14, 10, and 3 Gyr, followed by quiescent fading (see 
text for details). 
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Fig. 3. — The distribution of A/K vs Ha for all the galaxies in this sample. Typical errors in 
A/K are ~ 0.03 to 0.1. Shown in red are model tracks for the 3 models shown in Figure 2. 
The K+A galaxy population was selected with the cuts shown in green. 
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Fig. 4. — The distributions of the ratio of the two star formation indicators, Ha EW and 
A/K, for four subsamples of samplel2, selected to have A/K > 0.2, 0.3, and 0.4. Note the 
evidence for the existence of a separate population of K+A galaxies around zero. 
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Fig. 5. — Spectra of members of a randomly selected subsample of the K+A population, 
selected from A/K bins so as to span the range. The grey lines show the data and the black 
overplotted lines show the "K+A" best fit. 




Fig. 6. — The luminosity function — i.e., luminosity distribution of the sample weighted by 
1/Vmax — for the entire sample and the K+A galaxies. The plotted uncertainties are simply 
the l/\4iax weights added in quadrature (divided by the bin width). 




Fig. 7. — The distribution of A/K, weighted by l/Knax, to make a "ratio function" . Overplot- 
ted are curves proportional to |dt/d(A/K)| for the the "10" and "3" Gyr models of Figure 2. 
These curves are equivalent to "steady-state" predictions made under the assumptions that 
(a) all galaxies have the same star-formation history, and (b) galaxies are passing through 
the A+K phase at a rate that varies slowly on Gyr timescales. The plotted uncertainties are 
simply the 1/Knax weights added in quadrature (divided by the bin width). 
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Fig. 8. — The top left panel shows the distribution of luminosity density in galaxies in the 
sample in the color-Sersic-index plane. The color is measured in blueshifted SDSS band- 
passes ^'^g and ^'^r, and the Scrsic index n is found by fitting to the radial profile in the ^'^i 
band. Each galaxy data point has been weighted by the ratio of its luminosity Lo.ij in the 
band to its comoving selection volume Knax, so the distribution of weight in the plot is pro- 
portional to comoving luminosity density. The contours enclose 52.0, 84.3, and 96.6 percent 
of the total luminosity density. In this panel, galaxies separate into disk-dominated (blue, 
low n) and bulge-dominated (red, high n) populations, marked by crosses. The other panels 
show similar plots, but for different sub-samples, cut on A/K and HaEW. The "K-l-A" panel 
shows those galaxies making the high A/K cut and the low HaEW cut. Galaxies selected 
to have low A/K or low HaEW appear to be concentrated (bulge-dominated); galaxies se- 
lected to have high A/K or high HaEW appear to be exponential (disk-dominated). The 
K-f-A galaxies have the colors of disk-dominated galaxies but many have the radial profiles 
of bulge-dominated galaxies. 
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Fig. 9. — Similar to Figure 8 but showing the distribution of luminosity density in galaxies in 
the sample in the color-surface-brightness plane. For each galaxy, the surface brightness /io.ij 
is the mean surface brightness inside the Pctrosian R^q radius of the best-fitting Sersic radial 
profile in the ^'-^i band. In this panel, again, galaxies also separate into disk-dominated (blue, 
low ii) and bulge-dominated (red, high n) populations, marked by crosses. Galaxies selected 
to have low A/K or low HaEW appear to be high in surface brightness (bulge-dominated); 
galaxies selected to have high A/K or high HaEW appear to be lower in surface brightness 
(disk-dominated). The K-l-A galaxies have the colors of disk-dominated galaxies but many 
have central surface brightnesses as high as, or even higher than, bulge-dominated galaxies. 
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Table 1. average overdensities around galaxy subsamples 



subsample 




N 


>ii/jtotal 


(Si) 




all 




125255 


1.000 


19.8 ± 1.7 


1.86 ±0.19 


disk-dominated {n 


< 2.0) 


49615 


0.433 ± 0.007 


13.7 ± 1.0 


1.55 ±0.19 


bulge-dominated {n > 2.0) 


75640 


0.567 ±0.007 


24.5 ±2.2 


2.10 ±0.18 


K-l-A galaxies (see 


text) 


1143 


0.009 ±0.001 


15.0 ±2.3 


1.64 ±0.20 



^Fraction of the total galaxy luminosity density in the band (i.e., total of 
-^o ij/Knax) contributed by each subsample. 



